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ABSTRACT 

, We present a generalization and refinement of the Sachs- Wolfe technique which 

ON I unifies many of the approaches taken to date and clarifies both the physical and the 

mathematical character of the method. We illustrate the formalism with a calcula- 
c3 ■ tion of the behavior of light passing a moving lens on a Minkowski background. 



Subject headings: gravitation — cosmology: cosmic microwave background — cosmology: 
gravitational lensing 



j> ■ I. Introduction 

^ \ A perfectly homogeneous and isotropic cosmic microwave background radia- 

tion (CMBR) is possible only in a Friedmann-Lemaitre-Robertson- Walker (FLRW) 
spacetime (Ehlers, Geren and Sachs, 1968). For this reason, the extreme isotropy 
of the CMBR has often been taken to indicate that our Universe is, in some sense, 
well approximated by one of the FLRW models. Nevertheless, there is structure in 
■ the Universe that must induce distortions into the CMBR. 

In the usual scenarios, the structure that we observe today evolves via gravi- 
tational instability from small inhomogeneities of the matter fields at the epoch of 
matter-radiation equality, which occurs at temperature Teq ~ 5.5fi /i 2 e V (where 
h = -ffo/lOOkms -1 Mpc -1 ), when the formation of galaxies and clusters by gravi- 
tational collapse begins. (Peebles and Yu, 1970; Bond and Efstathiou, 1984; Bond 
and Efstathiou, 1987; Suto, Gouda, and Sugiyama, 1990). An important prediction 
of gravitational instability models was first noticed by Sachs and Wolfe (1967) who 
showed that even small amplitude perturbations of FLRW spacetimes at recombi- 
nation (temperature Trec ~ -31 eV, notably close to the epoch of matter-radiation 
equality) cause appreciable temperature fluctuations in the CMBR. 

There are two main problems associated with a Sachs- Wolfe calculation; relating 
the metric perturbations to interesting physical perturbations of the matter fields, 
and determining the behavior of null geodesies in the perturbed spacetime. Solving 
the first problem amounts to constructing the physical model to be investigated. 
The behavior of electromagnetic radiation in the model can then be understood 
in the geometric optics limit by examining the behavior of null geodesies (Jordan, 
Ehlers, and Sachs, 1961). The original work of Sachs and Wolfe (1967) concerned 
itself with a particular perturbed Einstein de-Sitter spacetime and with a single 
observable property of its null rays, the redshift. This is an important quantity 



because it is directly, and simply, related to the temperature profile of a thermal 
radiation field such as the CMBR. Much subsequent research has been devoted 
to constructing perturbations about Einstein de-Sitter spacetime which embody 
specific physical characteristics (e.g. spheroidal perturbations) and using the equiv- 
alent of the Sachs- Wolfe formula appropriate for the specific perturbations in order 
to understand the temperature pattern of the background radiation in these models 
(Grishchuk and Zel'dovich, 1978; Linder, 1988a; Argiieso and Martinez-Gonzalez, 
1989; Argiieso, Martinez-Gonzalez, and Sanz, 1989). A formula for the redshift of 
null geodesies valid for perturbations about the curved FLRW models as well as the 
flat was produced by Anile and Motta in 1976. 

Since gravitational lensing may be thought of as an aspect of the theory of 
perturbed null geodesies, some researchers have used Sachs- Wolfe like calculations to 
investigate the bending of light rays. Linder has considered the spatial components 
of the perturbed geodesic equation on an Einstein de-Sitter background, in addition 
to the timelike component, and has analyzed questions usually found in the realm of 
gravitational lens theory with a Sachs- Wolfe type calculation (Linder, 1988b, 1990). 
Martinez-Gonzalez, Sanz, and Silk (1990) have exploited the fact that, under certain 
assumptions, metric perturbations of the flat FLRW spacetime may be expressed 
by a small, scalar potential which obeys an expanding-space analog of the usual 
Poisson equation. They use this result to write formulae for both redshift and 
deflection angle in physically transparent forms. 

A complete generalization of the Sachs- Wolfe formalism should solve for the 
timelike and spacelike components of null geodesies in a metric perturbed spacetime 
of arbitrary background. Such a generalization is presented here for the first time. 
When the background is taken to be Einstein-de Sitter our technique agrees with 
that of Sachs and Wolfe for the expression of the photon redshift and with that of 
Linder (1988b) for the perturbed photon path. When the background is taken to be 
a curved FLRW spacetime, our results agree with those of Anile and Motta for the 
photon redshift (the quantity with which they were concerned), but our method, 
unlike theirs, can also be used to describe the spatial components of the perturbed 
photon wavevector or lensing. Their formalism is inadequate for these issues because 
the spatial components of their primary equation (equation (7) of Anile and Motta, 
1976), do not take into account the distinct nature of the background and perturbed 
paths to the correct order. 

Since our formalism applies to general backgrounds it is able to handle interest- 
ing cases which were not previously amenable to a Sachs- Wolfe type calculation. For 
example, our method may be used to investigate perturbations of a Schwarzschild 
background by long-wavelength gravitational waves. We also expect that current 
developments in mathematical cosmology will benefit from our more general for- 
mulation of the Sachs- Wolfe technique. For instance, to date work has focused on 
models where the density perturbations are small, in the sense that the fractional 
density fluctuations dp/p <C 1. However, recent work by Futamase (1989) and Ja- 
cobs, Linder, and Wagoner (1992) allows a description of a spacetime with large 
density contrasts by a small metric perturbation. While the background spacetime 
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is not necessarily FLRW, "clumpy" cosmologies described in this way could be excel- 
lent models of our Universe. In addition, the notion of an inhomogeneous spacetime 
behaving, "on average" , like an FLRW spacetime is not rigorously understood and 
there is no reason to believe that future models will behave like linear perturbations 
of Einstein-de Sitter space. We also hope that the technique presented here will be 
able to clarify a difficult point in lensing theory, the nature of the distance factors 
in the lens equation. Because of the difficulty in working with the rigorous optical 
scalar equations (Sachs, 1961), the relationship between the formal mathematics of 
lensing and the, very successful, use of the lens equation is not yet clear (Futamase 
and Sasaki, 1989; Watanabe and Tomita, 1990). The relationship of our technique 
to the optical scalar equations is the subject of a forthcoming paper (Pyne and 
Birkinshaw, work in progress). 

This paper is organized as follows. In section II we formulate an equation 
describing the local separation of geodesies of the perturbed metric from geodesies 
of the background. We then show that this equation is formally integrable as a 
power series in the curvature of the background spacetime. In section III we use 
our solution to find the behavior of light passing near a moving lens on a Minkowski 
background, obtaining results in agreement with those obtained to date using other 
approaches. In section IV we summarize our results. 

II. Null Geodesies in Perturbed Spacetimes 

(a) Conventions. In the following, Greek indices run over 0,1,2,3 while Latin 
indices run over 1,2,3, the zeroth component is time-like, and the summation con- 
vention is assumed. The metric signature is +2 and we take G = c = 1. The 
Riemann tensor convention is defined by equation (A4). 

(b) A Perturbed Jacobi Equation. Our starting point is the usual formalism 
of metric perturbation theory. We consider the true metric to be the sum of a 
background metric gffj, and a small perturbation, h^, 



9^,v 9jiu 




(1) 



We define g» v and g^ v by 
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It is then a simple matter to verify that to order h 
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g"" = g {0) ^ - (3) 



where h^ u is denned by 



h V = gW-OgQWhap . (4) 



As usual (in a metric perturbation theory), we raise and lower all tensor indices 
with the background metric. 

The Levi-Civita connection of 



r M a/3 = ^ MfJ (g*i3, a + g a *,t3 - g a /3,a) (5) 



can be split into zeroth and first order components in h 



r% = r% + r(% (6) 



where 



r (0)/i o _ 1 n (0) f ia f(0) , (0) (o) \ m 



and 
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The semicolon above and in what follows denotes covariant differentiation with 
respect to the Levi-Civita connection of gj$. We note that r^ 1 ^ is tensorial, at 
least in the usual restricted sense of metric perturbation theory (that is, under 
infinitesimal co-ordinate transformations). 

Let x^ fl (X) be a geodesic of the background spacetime with A affine. We will 
sometimes refer to x^ fI (X) as the "unperturbed path". It satisfies the geodesic 
equation in the unperturbed spacetime, (• = d/dX) 



x<°>" + r(°Kp (x {0) ) ± (0)a ± (0)/? = o. (9) 



The notation I^ ^^ (V 0) ) is meant to convey explicitly that the connection terms 
are evaluated on the unperturbed path. 

Now consider the expression 



x ^(X) = x^(X)+x {1 ^(X) (10) 



where x^°^(X) is the unperturbed geodesic. At this point, both x M (A) and x^ fl (X) 
are unspecified and so equation (10) can be considered to define either x^ or x^ 1 ^ 
once the other is given. We will derive conditions on x^ fl (X) which will be necessary 
and sufficient if x M (A) is to be geodesic in the perturbed spacetime. In the following 
derivation, we will truncate three expressions by discarding terms of quadratic order 
in the products of x^ 1 ^ and x^ 1 ^ . We will often refer to this as "working to first 
order." The consistency conditions for the truncations made will be discussed at 
the end of the derivation. 

Differentiating equation (10) twice and using equation (9) gives 



5 



£M = _ r (0)M a/3 ( x (0)) ± (0)« ± (0)/3 + ^d)M. 



(11) 



On the other hand, if x M (A) is to be an affinely parameterized geodesic of the 
perturbed spacetime we must have 



= _ r (°)M a/ , ( X ) (±(0)^(0^ + 2±(°) a ±( 1 ^) (12) 



where we have used equations (6) and (10), keeping terms only to first order. Pro- 
vided no singularities intervene (which we can ensure by keeping to small enough 
neighborhoods and regular points), the connection terms near the perturbed path, 
x, may be expanded about their values on the unperturbed path, x^°\ as 



(x) = r< >" 0/9 (*«») + r^ a/3 , r (*<°>) + 
r%(x) = r%(/>) + ... . 



(13) 



Substituting these expansions into equation (12) gives 



(0)0 



(14) 



Comparing equations (14) and (11) we conclude that x M (A), defined by equation 
(10), will be an affinely parametrized geodesic of the perturbed spacetime provided 
x^ 1 ^(\) satisfies the system of four coupled, second-order differential equations 



(^+^H* (1,=/ ' (15) 
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where the 4x4 matrices A and B and the four-vector / are defined by 



Wa = r (0)H Taak (0)r k (0)a (16) 

jh _ _r (1)At T(J A; (0)r /c (0)fJ 



Anticipating the interpretation of null geodesies as photon paths we have written 
as fc(°). A 4 x 4 identity matrix is the implied coefficient of the second-order 
derivative in equation (15). We are using the matrix notation as a convenient 

shorthand for unambiguous summations, e.g ^Bx^ 1 ^ = B^i^ 1 ^. We emphasize 
that equation (15) is to hold along some segment of the unperturbed path and 
generates solutions for the separation, x^\ of the perturbed path relative to the 
background path. 

We still need to discuss where the above equation is actually valid, that is, the 
appropriate consistency criteria for the approximations made. We will use Futa- 
mase's (1988) goodness-of-fit parameters for this purpose. Let the typical magnitude 
of a component of h relative to that of a component of </°) be written e 2 . Write 
the typical scales of h and as / and L, respectively. Put k = l/L. Also, for 
simplicity (though without loss of generality), assume that the initial conditions 
appropriate for the sought-after solution are given at A = 0. Let k^ denote the 
magnitude of a typical component of the unperturbed wavevector at A = 0. The 
consistency criteria for the approximations made above may then be found with the 
help of equations (9), above, and (20) below to be e 2 1, e 2 k, and 



e 2 Llnl^ + lj «min(/,L). (17) 



This last condition may be replaced with the stricter e 2 Xk^ <ti min(/, L). 

In appendix A we show that equation (15) is simply a perturbation of the Jacobi 
equation of the background spacetime. In its most elegant form, then, our equation 
reads 



Vl (0) x^-R(k (0 \xW)k^=f (18) 
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where V is the Christoffel connection of the background spacetime and R is its 
curvature tensor (with implicit superscript (0)). Remembering the tensorial nature 
of /, the tensorial character of our equation is now immediate. The physical meaning 
of the equation is also clear. The geodesic of the perturbed metric differs from that 
of the background because of 

(1) a force-like term associated with the perturbation (encoded in /) and 

(2) the usual geodesic deviation induced by the background spacetime. 

(c) Comment on Initial Data. We now address the issue of what are appro- 
priate initial data for the perturbed Jacobi equation. We know that the geodesic 
which we seek to construct is uniquely specified by a point of that geodesic, p, and 
the tangent vector to the geodesic at that point, k p . As equation (10) indicates, we 
can construct a segment of the desired geodesic locally about p from any geodesic 
of the background which (in the naive sense) passes close to p with tangent vector 
close to k p . An important point, however, is that the initial data specifying the 
desired geodesic does not, in general, generate a null geodesic of the background 
because k p is generally not null in both the perturbed and background metrics. As 
we wish to construct null geodesies of the full metric using a null geodesic of the 
background, we conclude that the initial condition k^(X p ) = is not appropriate 
for our equation. However, in the regions far from a localized perturbation (where, 
loosely, h — > 0) this subtlety is entirely avoided. 

It is instructive to examine this point in another way. Suppose that we have cho- 
sen a null geodesic of the background and that a solution, x^^A), of the perturbed 
Jacobi equation has been found. The condition that the constructed geodesic be 
null, to order h, in the perturbed spacetime is 



g^(x)k»k u = h^(x^)k(°^k {0) » + 2g$(x {0) )k {0)fl k^ 1/ 

+ 2g^Jx^)x^k^k^ (19) 

= 0. 



Making use of the co-ordinate invariance of scalar quantities, we evaluate the last 
term on the RHS of the first equality by choosing co-ordinates adapted to at 
any given point so that we can replace the ordinary derivative with the covariant 
derivative of g(°\ Metric compatibility then tells us that this term vanishes. We 
are left with the condition that everywhere along the background geodesic we must 
have 
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(20) 



We now show that this condition can always be enforced. We choose to work 
with co-ordinates for which the background connection coefficients vanish along the 
unperturbed path. In these co-ordinates, equation (15) implies 



(21) 

where the first term on the RHS vanishes as a consequence of the identity 

^L^^V^ = o 



written out in the choosen co-ordinates. The RHS of equation (21) may be shown, 
with the aid of equation (8) , to be 



= _i fe (oH(o)- fc (o)/j ( ha ^ _ r (°)%v - r<®%f,h ap ) . (22) 



We have thus shown that 



d_ (,(0),(0H(iV) = R. (-\ Kv m^) (23) 



where D/dX is the covariant derivative of g(°> along x^. 

On scalars we have D/dX = d/dX so we can conclude that if equation (20) is 
satisfied at any point along the background geodesic it will be satisfied at every 
point. Notice that it is very simple to impose equation (20) at the initial point but 
that, in general, k^ 1 ' = is not an acceptable solution because at a given point a null 
vector in the perturbed metric is not necessarily null in the background metric. The 
calculation above thus amounts to a direct check that (k^ + k^ l \ k^ + k^') g (o) +h 
is a constant of geodesic motion to order h, as it must be if this procedure is to be 
meaningful (the angle brackets denote the inner product with respect to the metric 
of the subscript) . 



9 



(d) A Solution to the Perturbed Jacobi Equation. Having found the 
domain for which our equation, (15), is valid and having shown that the equation 
generates null paths from null paths, we construct a formal solution. The techniques 
we will use are simply the matrix analogues of familiar methods for dealing with 
ordinary differential equations. First we perform a change of variables in order to 
eliminate the first derivative term in our equation. To accomplish this, let P(X, a) 
be a 4 x 4 matrix function of two real arguments, and v(X) a vector such that 
x^ = Pv. Then in terms of v, equation (15) becomes 



v + P' 1 (2P + AP) v + P' 1 (P + AP + BP) v = P" 1 / (24) 



where we have assumed P non-singular (justified below). Now choose P to satisfy 



p =-^ p > ( 25 ) 



The solution is the path-ordered exponential, 



P(X, a)=V exp (~ £ A{t) dr\ (26) 



(a short introduction to the path-ordering symbol, V, is provided in appendix B). 
We remark that P is exactly Synge's parallel propagator (Synge, 1960). We will 
call P the connector, following De Felice and Clarke (1990). Anticipating our final 
result, we see that the effect of the change of variables above is basically to untangle 
the co-ordinate basis, just as if we were working with a parallel propagated tetrad 
frame. 

In writing equation (26) we have set P(a, a) = Id, the four-dimensional identity 
matrix. This initial condition and equation (25) can be used to show the important 
identity 



P(X,X 1 )P(X 1 ,a) = P(X,a) 
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(27) 



and so 



P(X,a) 



-l 



P(a, A) 



(28) 



proving that P is invertible. 

With P chosen in this manner, equation (24) now becomes 



v + P 



-l 



( 




(29) 



Written out explicitly using equation (16) the quantity in parentheses reveals itself 
to be -R^upak^k^f, which we will write as the matrix —R, using notation 
consistent with (18). 

We write (29) as a first-order system 



which can be solved by constructing the transition matrix of the system, U(X, a), 
which is the solution of its homogeneous part (Humi and Miller, 1988). This is the 
familiar method of Greens function solution. We will use the terms transition matrix 
and Greens function interchangeably in refering to U(X, a). Another path-ordering 
is needed to actually construct the transition matrix, with the result 



where the boundary condition U (a, a) = Id has been imposed. As a result, equations 
(27) and (28) remain valid with U written in place of P. Defining the column vectors 
y = (v,v) and s = (0, P~ l f) the formal solution to equation (29) with initial value 
boundary data is 




(30) 




(31) 
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y(X) = U(X, a)y(a) + U(X, r)s(r) dr 

J a 



(32) 



Confirmation that this expression solves equation (29) is provided by differentiation. 

While (32) is a formal solution to the perturbed Jacobi equation (29; and hence 
15), at this point its utility is far from evident. In particular, the presence of two 
path-ordered exponentials makes the calculations difficult. Had we not made the 
initial change of variables from x^ 1 ' to v we could have solved the system without 
recourse to the connector. We would then need only a single path-ordering, that 
responsible for constructing the appropriate Greens function. Appendix C outlines 
this approach. We feel, however, that the strength of the particular transcription 
of the solution to the perturbed Jacobi equation given above, equation (32), lies in 
its explicit geometric character. This is most evident by choosing co-ordinates for 
which the Christoffel coefficients of the background vanish along the background 
geodesic. In this case equation (32) applies with P = la, s = (0, /), y = (x^\ k^ 1 ') 
and 



U(\,a) = VexpU (° R o d ) dr) (33) 



clearly expressing the relevance of the background curvature to the solution. 

III. The Moving Lens 

As an illustration of our method, we calculate the asymptotic behavior of a light 
ray passing by a point mass moving on a Minkowski background. Elements of this 
scenario have been analyzed by Birkinshaw and Gull (1983) without a fully general- 
relativistic treatment. Our formalism provides an easy and unified way of handling 
the problem in the context of general relativity. This particular case is especially 
simple because Minkowski space is flat, allowing us to choose co-ordinates such that 
A = B = 0, and hence allowing equation (15) to be solved by direct integration. 
In appendix D we show how the propagator techniques of section II can be used to 
gain an equivalent solution. 

(a) The Metric. Let primed co-ordinates denote the frame in which the lens 
is at rest at the origin of co-ordinates. In this frame the metric is of the usual 
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Schwarzschild form, 



ds >2 = _ ( 1 _ ^ d1 p + (l _ ^ 1 rfr '2 + /2 ^/2 + r /2 gin ^,2 ^ (34) 



where m is the mass of the lens. Expanding to first order in m/r' we break the metric 
up into a sum of the Minkowski metric in spherical co-ordinates and a perturbation 
of the form 

W = ^ (OS' + £4') • (35) 



Transforming to Cartesian co-ordinates puts the background Minkowski metric into 
its usual form, diag(— 1, 1, 1, 1), and puts the perturbation into the form 



W {x) = — 

ho'i' (x) = hi>o> (x) = (36) 
m 

Ji-Xi'Xji 



where 



r '= L'2 + ?/ /2 +/ 2l 1 /2_ ^ 



Note that the x^ appearing in equation (36) are not vectors but the primed co- 
ordinate functions. In particular no index raising or lowering via a metric is taking 
place, Xi' = x % . This remark will apply equally well to the unprimed co-ordinate 
functions. 

We are interested in the background metric and metric perturbation in the 
frame in which the lens appears to be moving with velocity v, with the underbar 
denoting a three-vector quantity. We will use unprimed variables for this frame. If 
we agree to work only to linear order in v the connection between our two sets of 
co-ordinates is simply a linearized Lorentz transformation 



i : =i - £ ; £ (38) 

X = x — vt. 
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where the dot product is a convenient shorthand for summation of spatial indices, 
e.g. v ■ x = v l Xi. Under this transformation the background Minkowski metric is 
unchanged while the metric perturbation becomes 



h 



2m 



oo 



2m 2m 
hoi = h i0 = Vi o-{v-x)Xi 

2m 

i^ij — \X{Xj XjV%t XjVjtj 



(39) 



with 



r = \J (x — vt) ■ (x — vt) . 



(40) 



The global Galilean transformation, (38), may be used to define the three- 
vector v because the background is Minkowski and because products 0(vh) are 
ignored, so that v l = Vi. Although this is not globally valid in general relativity, the 
approximations adopted are consistent to the order claimed, as may be discovered by 
repeating the calculation using the weak-field equations of general relativity, where 
the perturbation can be written diag(— 20, —20, —20, —20) with the Newtonian 
potential of the perturbation, taken to be a moving point mass (see e.g. Weinberg, 
1972). 

(b) Connection Terms. Simple calculation using equations (8) and (39) 
produces (to 0(v,h)), 



r (i)o _ 
1 oo — 


m 


x 


r (i)o 

1 Oi _ 


m . 


- Vit) 


r (l)0 _ 
L jk ~ 


2m . 


' ^> °i k ~ ^3 (¥* + v kXj) --^-(v-x) XjX k 


r (l)i _ 

1 oo — 




-vH) 


r (l)i _ 

L 0j - 


3m . 

7^ 


s m , 

X ) X'lXj ~\~ 2 {^XjX)^ C>Xj%)jj 


L jk - 


2m 

73" ( Xi 


3m 

f^it) Oj'fc g - yXiXjXk XiX jVj$t XjXfcVit XiXfcVjt) 



(41) 

(c) The Unperturbed Path. Let t = be the time of closest approach of the 
photon and lens, and use co-ordinates centered on the lens with co-ordinate axes 
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chosen so that at t = the photon wavevector, n , points in the y direction and the 
photon position vector, r , lies along the z axis. With this choice of co-ordinates 
the unperturbed path of the photon is simply = (A, 0, A, ro) where ro = |r | is 
a constant of the motion. 

Analysis of the consistency criteria for our method quickly reveals that we may 
work over the entire unperturbed path provided e 2 < 1 and e 2 </t. The first of 
these inequalities limits us to regions of spacetime where the Newtonian potential 
of the perturbing mass is small, and the second restricts the mass to move at non- 
relativistic velocities. The second of these conditions we have already imposed. The 
first amounts to restricting our attention to weak lensing scenarios, that is, lensing 
for which the impact parameter is much larger than twice the Schwarzschild radius 
of the lens. 

(d) The Bend Angle. We will call the deflection angle in the yz-plane the 
bend angle. In the small angle approximation it is given by 



_dz 

dy 

z 

y 



,(D ( 42 ) 

IT? 1 ) 

J -co x ' 



(using equation (15) in its first integral form, with A = B = 0). Substituting the 
co-ordinate values on the unperturbed path into the expressions for the connection 
coefficients, switching variables to s = A/ro, and working to first order in v gives 



m f+°° 

Pbend = / 

r J-oo 



ds 



-oo (l + S 2)5/2 



, ^ q 15s (v z + SVy) 

3(1-2^) -v z s + 2v z s 3 + v y> 



1 + s< 



(43) 



which may be integrated trivially to yield the result 



4.777, 

#bend = (1-V y ) (44) 

7*0 



or, in three-vector form, 



4tm 

#bend = (1-y-no). (45) 

^0 
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which is valid to first order in m/r and v. This is the first result of our method 
and a short calculation shows that it agrees with a Lorentz transformation of the 
usual (static) lens deflection angle to a frame in which the lens moves with velocity 
v. 

(e) The Frequency Shift. The frequency shift of the photon between emis- 
sion, e, and observation, o, is defined by 



v v e (46) 

= (l + ^)- 1 -l. 



where the redshift, z, is given by 

(k ■ u) 



(1 + *) _1 = 



(k ■ u) e ' 



with u and u e the four- velocities of the observer and emitter, respectively. Keeping 
terms only to linear order, 



-i _ k^\o)u^{o) + k^(o)u^(o) + k^\o)u^(o) 
~~ *i 0) (e)it(°)M(e) + 4 0) (e)w( 1 )M( e ) + k^ (e)u(°)<*(e) 



We take the observer and emitter to be far from the lens and at rest relative to each 
other (letting the observer and the emitter have some relative motion would merely 
give rise to the usual Doppler terms). To be precise, the observer and emitter four- 
velocities are parallel translates of each other along a path passing far from the lens, 
this relationship being, at least asymptotically, path-independent. This allows us 
to write u^°\o) = u^°'(e) = (1, 0, 0, 0). Furthermore, in the asymptotic limit we are 
considering, u^^o) = u^^e) = (0,0,0,0). Making these substitutions in equation 
(47) and combining the resulting expression with equation (46) yields 



Ai/ _ki°\o) + k^(o) L (48) 



^ k^(e) + k^(e) 



We showed in section II that for a perturbation of restricted scale we could self- 
consistently assume that the photon is unperturbed at emission, so that fcW(e) = 0. 
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In addition, = g^k^ ^ = —i^ = — 1 at both emission and reception, from 
which we deduce 

— - -k {1) (o) 
v - k ° [0) (49) 

= k^°(o)=i {1) (o). 



Taking our emission point to be A = — oo and our reception point to be A = +00 
equations (49) and (15) lead to 



— = - r ( r o»" + 2r s'" + r -°) dx < 5 °) 

V J — 00 



Substituting in the co-ordinate values for the unperturbed path, transforming vari- 
ables to s = A/ro, and performing a few algebraic manipulations leads to 



Au m f+°° ds r „ , „ 9 n a ^ . 

= / rjK v z + s(2- 3v y ) - 6s 2 v z + 2s 3 1 - 5v y ) 

v r Q J-oo[i-2sv z + (l-2v y )s 2 ]° /21 

from which an expansion to 0(v) yields 



(51) 



Au Am 

= v z 52 

v r 



or, in three-vector notation 



— - -^v-r Q . (53) 



v 







Like the formula for the bend angle, this result agrees with previous calculations of 
the frequency shift (Birkinshaw and Gull, 1983): to 0(v) the frequency shift arises 
from transverse motion of the lens. 

(f) The Skew Angle. We will call the angular deflection out of the yz-plane 
the skew angle. The special relativistic treatment of the scattering of a photon by a 
moving point mass represents an instantaneous interaction between the photon and 
the lens and so, necessarily, gives a vanishing skew angle. Our method accounts for 
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the interaction of the photon and the lens over the entire photon path and confirms 
that the skew angle vanishes by direct calculation to order mv/ro. 

In the small angle approximation the skew angle is given by 



ft - — 

"skew , 

dy 

x 

y 



X 



(1) 



1 + 

(1) - 



X 



(54) 



r + OO 

= - r 

J — oo v 



(1)1 , 9F (1)1 
00 ' ZL 02 



+ 4f 



)d\. 



which becomes 



^skew 



+oo 



s ds 



r () J-oo [l _ 2v z s + (1 - 2v y ) s 2 f /2 
3mvT r + °° s 3 ds 



ro 



(55) 



oo [l_2z; 2 s + (l-2^)s 2 ] 5/2 ' 



which vanishes to first order in v. 



IV. Summary 

We have described a method for constructing null geodesies in arbitrary met- 
ric perturbed spacetimes using null geodesies of the background spacetime. The 
method fully generalizes the usual Sachs- Wolfe technique for calculating tempera- 
ture fluctuations of the CMBR in metric perturbed spacetimes. Because our method 
constructs both the spatial and timelike components of the perturbed geodesic it is 
able to address questions of interest in gravitational lens theory, such as the bend 
angles of the true path relative to the unperturbed path. We have provided an 
explicit illustration with the calculation of the behavior of a photon passing a mov- 
ing lens. A forthcoming paper will show how our method can be used to calculate 
other quantities of gravitational lens theory, such as the amplification undergone by 
a bundle of light rays passing a given perturbation. 



Appendix A: Rewriting the Jacobi Operator 

We prove the equivalence of equations (15) and (18). This short calculation has 
an interesting history. It is implicit in Weinberg (1972), section 6.10. We produced 
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the following proof in the course of trying to gain a physical understanding of 
equation (15). After we had done so, a comment in Burke (1985) drew our attention 
to a paper by Faulkner and Flannery (1978) where the same calculation appears in a 
different context. It appears that this calculation is either obvious or not depending 
on who is presenting it. It was not obvious to us, and we feel its importance to the 
physical understanding of equation (15) warrants its inclusion here. 

Let D/dX denote covariant differentiation along the curve with the back- 
ground connection. Then for an arbitrary vector v 



and 



Using the geodesic equation for this becomes 



^%M = ^ + 2r (0W0),^ 



dX 2 ^ dX 2 af3 dX 

+ T^k^k^v? + Y^Y^k^k^vP ( A3 ) 

- Y^r^ a k Wa k Wp v p 



The Riemann tensor of the background is given by 



d(0) M ( r (0)i* _ r (o)/x r (o)M r (o) P _ r (0) Mr (o)/A (M \ 



so that (A3) can be written 
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Comparing this equation to equation (15) completes the proof. 



Appendix B: The Path-Ordering 

This appendix provides a brief introduction to the use of the path-ordering 
symbol. Recall that in section II we needed to solve a matrix system of the form 



P = MP (SI) 



subject to P(a, a) = 1^. If the quantities P and M were functions, the solution of 
this system would be the usual exponential. An exponential with matrix argument 
will not work, however, for the reason that a matrix function and its derivative 
matrix do not, in general, commute. 

Rewriting equation (Bl) as the integral system 



P(X, a) = l d + I* M(r)P(r, a) dr (B2) 

J a 



allows a possible solution, order by order in M, to be written down by iteration, 



P(X, a) = l d + ( X M(r) dr + ^ M(r) dr f M(r') dr' + . . . (S3) 

J a J a J a 



A simple proof of convergence then shows that this expression is, in fact, a valid 
solution. The region of integration for the n-th order term (in M) is known as the 
n-simplex (the n-dimensional analog of the triangle). We can extend the region 
of integration to the n-cube in two steps. At each order we must symmetrize the 
n-fold product M(t)M(t') . . . M(r' - ") completely in its arguments. This extends 
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the argument of the n-simplex integral to an argument defined over the entire n- 
cube and ensures that each of the n! simplices in the n-cube contributes the same 
amount to the total integral over the n-cube, an amount equal to the value of the 
original integral over the single simplex. Our second step, then, is to divide each 
term by its overcounting factor, n! at n-th order. These combinatorial factors yield 
the exponential. Thus the path-ordering symbol amounts to a notice to perform an 
integration written over a cube only over its "lowest" (in the sense of the discussion 
above) simplex. 

Appendix C: Another Transcription of the Solution 

In this appendix we present another form for the solution of the perturbed 
Jacobi equation, needing only a single path-ordering but containing terms whose 
geometrical meanings are more obscure than those in equation (32). We first write 
equation (15) as an eight-dimensional first order system, 



where Id is a 4 x 4 identity matrix and A, B and / are given by equation (16). 
From this point we will denote the 8x8 matrix which is the coefficient matrix of 
the associated homogeneous system by M. 

We can now proceed to solve equation (CI) in exactly the same manner as we 
solved equation (30) in section II. We first obtain the associated transition matrix, 
U, which solves 



subject to U(a, a) = Id- We note that the equation for U is solved by a path-ordered 
exponential 





d_ 

dX 



U(X,a) = M(X)U(X,a) 



(C2) 




(C3) 



Having obtained U it is easy to check by straightforward differentiation that 
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equation (CI) is solved by, with y := (x^\ k^) and s := (0, /), 



y(X) = U(X,a)y(a) + / U(X, r)s(r) dr 



(C4) 



This is the solution we desired to obtain. 

The equivalence of this solution and that given in the main body of the paper 
is easily established by working in the co-ordinate system for which the Christof- 
fel coefficients vanish along the unperturbed path. The different appearance of the 
two solutions is essentially equivalent to the two different ways of writing the Jacobi 
equation. The usual expression for the Jacobi operator, involving covariant deriva- 
tives and the Riemann tensor, has the advantage of being an obviously geometrical 
quantity. Writing the Jacobi operator as the LHS of equation (15) obscures its 
geometrical meaning but requires fewer computations of coefficients. The solutions 
mirror the strengths and disadvantages of the two starting formulations although 
they are completely equivalent. 

Appendix D: The Greens Function of Minkowski Space 

In this appendix we construct the Greens function appropriate to a perturbed 
Minkowski space and show its equivalence to direct integration of equation (15). In 
co-ordinates for which the connection terms of the Minkowski background vanish, 
the connector is the identity and the transition matrix can be calculated instantly 
from equation (33) with R = 0, yielding 



Simple matrix multiplication verifies that this solution satisfies conditions (27) 
and (28) (rewritten with U in place of P). It is easy to see that the bottom four 
rows yield a solution for the photon wavevector equivalent to simple integration 
of the RHS of equation (15). To see that the photon paths calculated by the two 
techniques agree, start from equation (32). Assuming vanishing initial data for 
convenience and using the form of the propagator given in equation (Dl) yields 




(Dl) 





22 



On the other hand, a direct second integration of equation (15) gives 



x^(X)= f k^\r)dr=f dr f s{r')dr' (D3) 

J a J a J a 



which can be seen to be equivalent to (D2) after integration by parts. 



We would like to thank Sean Carroll for many helpful discussions, and the 
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